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Abst ract 



In this report we list the Fortr^xn subroutines for 
evaluating the confluent hy pe rgeomet r i c functions M(rt,b;x) and 
U(a,b;x). Ifiese subroutines use the stable recurrence 
re 1 at ions given e . g . in Wimp. 
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I ntroduct: ion 



It: is well known t^hat the ordinary d i -f -f e rent i al equation 






+ ^ - ay = 0 



dx2 dx 



las a solution 



y(x) = AM(a 5 l;x) + BU(a,l;x) 

it a is not a negative integer. 

This problem arises e.g. when solving the linearized 
shallow water equations with the Tull linear variation in depth 
included (see Williams, StaniTorth and Neta, [1]) • 

The computation oT the conTluent hypergeometric Tunctions 
is based on the Miller algorithm (see e.g. Wimp, [2]). In 
general, one has a second order diTTerence equation 

z(n) + a(n)z(n + l) + b(n)z(n+2) =0, n > 0, t>(n) 0 . 

IT b(n) = 0 Tor some n. in some cases one can make a change oT 
variable V(n) = A(n)z(n) which will produce an equation oT the 
desired type. Let w(n) be a nontrix^ial solution and the sum oT 
the n ^rmalizing series 



oo 



s = E c(k)w(k) ^ 0 
k=0 
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is known. Let N be a large integei' and define sC|^(n) . 
N+1 . by 



Zn C'^) 



rO n = N+1 

ll n = N 



Z|M(n) + a(n)z^(n + l) + b ( n ) (n + 2) =0. n = N-1. . 



One can approximat.e w(n) b^' W|,g ( n ) 



w^(n) = Szf^(n)/.S|^ 



where 



N 

Sn = E 

k=0 



c(k)Zf^(k) . 



The algorithm is said to converge iT 



'' N ( — CC . 



Ihe function M(a.b:x) satisfies the recui rence i- 



(2n+b+2) (n+a)z(n) 



(2n+b) (2n+b+2 

X 



- ( 2n + b ) ( n-t-b-f-l -a) z ( n-{-2 ) = 0 . 



The minimal solution is 



w( n ) 



x^a)n 

2n 



.M ( a-i-n . 2n4-b : X ) 



0 < n < 



... 1 . 0 . 



e 1 a t i o 1 1 
-|z ( n-f 1 ' 
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where 



(c)n 



r (n+c) 

r(c) 



The normalization relationship used in our subroutine is 

oo ( 

S = b-1 = g ^ , / (b^l),^(b+2k-l)w(k) . 

k=0 

An obvious modification must be made if b = 1. The algorithm 
is not defined if b. b+l-a, a are negative integers or zero. 
The function U(a,b;x) satisfies the relationship 



(n+a) (n+a+1 -b)z(n) - (n+l)[2 (n+a+1 ) Tx-b] z (n+1 ) 

+ (n+l) (n+2)z(n+2) = 0 



The minimal solution is 



w(n) 



^ ( ®-) n ( ^+1 “ t>) n 

n 



U (a+n . b ; x) 



for I arg x| < tt . A normalization relation is 

1 = £ w(k) . 

k=0 

In tlie next section we give a listing of the Fortran 
subrout i nes . 
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Subroutine Miller 



SUBF^OUTI NE M I LLEFi (N . AF.PFFA . F3ETA . X , S . SS . COEFF) 
INTEGEFF N 

REAL*8 AURFA .BETA ,X ,SS 
REAL*S S(0:1000) 

EXTERNA F. COEFF 

C USES TFFE .J . C . F’ . MILLER ALGOR ITFIM TO COMF>UTE 
C S(0:N). 

C BEG I N 

INTEGER NN.F< 

REAL>t=8 T.D.EPS.A,B.C 
REALMS 0F.DS(0 : 1000) 

EPS = 0.000000001 
C INITIALIZE OLDS. 

DO 1 F< = 0. 1000 
OLDS(F<) = 0 

1 CONTINUE 

C CFFOOSE INITIAL NN . 

NN = N + 2 

C INITIALIZE F\. S AND T. 

2 F< = NN 

S(F<+1) = 0 
S(F<) =1 

CALL COEFF (F< . ALPHA . BETA . X . A . B . C) 

T = 2*C=+=S(F<) 

C TAF\'E A BACF<WARD RECURRENCE STEP AND UPDATE IT. 

3 F< = F< - 1 

CALL COEFF ( F< . ALPHA . BETA . X . A . B . C) 

S(F<) = A-S(F<+1) + B-S(F<+2) 

C CHECF< FOR OVERFLOW AND RESCALE IF NECESSARY'. 

D= DABS (S(F\) ) 

IF (D .GT. 1.D30) TIFEN 

C BEG I N 

CALL SCALE (K.NN.S.T.D) 

END IF 

IF (F\ .GT. 0) HHEN 
C BEG I N 

T = T + 2-C-S(F<) 

GO TO 3 
END I F 

T = 1 + C-S(O) 

DO 4 F< = 0. N 
S(F<) = S(F^)/^1 

4 CONTINUE 
TEMPORARY PRINT STATEMENT. 

PRIN1-. S(0) 

TEST FOR CONVERGENC'E. 

D = 0 

DO 5 F\ = 0. N 

D = D + S(Fx)»-2 
CONTINUE 
D = DSQiri(D) 

T = 0 



.5 
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C 

C 



999 



DO 6 K = 0, N 

T = T + (S(K) - 0LDS(K))**2 
CONTINUE 
T = DSQRT(T) 

TAKE ANOTHER STEP IF NO CONVERGENCE. 

IF (T .GT. EPS*D) THEN 
BEG I N 

NN = 2*NN 
DO 7 K = 0, N 

OLDS(K) = S(K) 

CONTINUE 

IF(NN .LE. 1000) GO TO 2 
PRINT 999, NN, ALPHA, BETA, X,T 
FORMAT(' ** NO CONVERGENCE ** NN AP CP X T 
END IF 
SS=S(0) 

RETURN 

END 



', I5,4E14. 
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SUBROUTINE COEFF (N . ALPHA . BETA . X . A . B . C) 

INTEGER N 

REAL*8 ALPHA . BETA . X . A . B . C 

COMPUTES COEFFICIENTS USED BY J.C.P. MILLER ALGORITHM FOR 
A CONFLUENT HYPERGEOMETRIC FUNCTION M(a.b:x) 

SEE JET WIMP. COMPUTATION WITH RECURRENCE REI.ATIONS. 

PITMAN 1984 I^P . 61 -C2 

BEGIN 

INTEGER M.K 
REAL*8.T.U. V,W 
S = 2*ALPHA - BETA 
T = N -1 ALPHA 
M = 2*N 
U = M + BETA 
V = U + 1 
W = V -I- 1 

A = (S/W + U/X)*V/T 
B = (N -l BETA - ALPHA + 1)*U/T/W 
J' = 1 

IF (N .GT. 0) THEN 
BEGIN 

S = BETA - 1 
DO 1 K = 1. N-1 
T = -T*(l+S/K) 

CONTINUE 
T = -T-(l+S/M) 

END IF 
C = T 
RETURN 
END 



SUBROUT I NE SCALE ( K . N . S . T . D) 
INTEGER N.K 
REAL*8 T.D 
REAL*8 S (0:1000) 

BEG I N 

INTEGER J 

1 = T/D 

DO 1 .1 = l\ . N 

S(J) = S(J)/D 
CONT I NT'i: 

RETURN 



END 



SUBROUTINE COEFU (N , ALPHA , BETA , X , A , B , C) 

INTEGER N 

REALMS ALPHA. BETA,X,A,B.C 

C COMPUTES COEFFICIENTS USED BY J.C.P. MILLER ALGORITHM FOR 

C A CONFLUENT HYPERGEOMETRIC FUNCTION U(a.b;x) 

C SEE JET WIMP, COMPUTATION WITH RECURRENCE RELATIONS, 

C PITMAN 1984 PP . 63-64 

C BEGIN 

INTEGER M.K 

REAL*8 S,T,U,V,W 

S = ALPHA + QFLOAT(N) 

T = S + 1 .DO 
U = S*(T - BETA) 

V = QFLOAT(N + 1) 

W = \' + 1 . DO 

A = (2=+=T + X - BETA)=+=V/U 
B = - V*W/U 
C = 1 
RETURN 

END 



Remark: The program "that: calls Miller must supply as a last 

parameter either COEFF (for M) or COEFU (for U) . 



8 



The subroutines are available on a diskette from eith^^r authoi 
upon request. These subroutines were tested ext^^nsively foi' 
vai’ious values of a. b and x. 

Remark : IT the parameter is a negative integc^r, the solution 

oT The differential equation is 



where n = -a. 

Ln(x) are Laguerre polynomials whose coefficients a- 
sat i sf y 



y = ALn(x) + B{ln|x|Ln(x) + g /^mX^} 



m :=0 



i 




1 



n 



- n 



The coefficients satisfy' 





m 



1 



n-1 



(m+l)2 




m 



n 






m - - 1 j 

2 

m 



m 



n-f I . n-f 2 . . . . 
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